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Abstract. We present evidence to suggest that the study of one dimensional 
quasi-exactly solvable (QES) models in quantum mechanics should be extended 
beyond the usual sl(2) approach. The motivation is twofold: We first show 
that certain quasi-exactly solvable potentials constructed with the sl(2) Lie 
algebraic method allow for a new larger portion of the spectrum to be obtained 
algebraically. This is done via another algebraization in which the algebraic 
hamiltonian cannot be expressed as a polynomial in the generators of sl(2). We 
then show an example of a new quasi-exactly solvable potential which cannot 
be obtained within the Lie-algebraic approach. 



1. Introduction 

Lie-algebraic and Lie group theoretic methods have played a significant role in 
finding exact solutions to the Schrodinger equation in quantum mechanics. In the 
classical applications, the Lie group appears as a symmetry group of the Hamilton- 
ian operator, and the associated representation theory provides an algebraic means 
for computing the spectrum. Of particular importance are the exactly solvable 
problems, such as the harmonic oscillator or the hydrogen atom, whose point spec- 
trum can be completely determined using purely algebraic methods. The concept 
of a spectrum generating algebra dates back to a paper by Goshen and Lipkin [1] , 
and was later rediscovered in the context of high energy physics by two different 
groups [2,3]. The study of spectrum generating algebras received much impetus in 
the subsequent years (see the review [4]) and it was soon applied also in the field of 
molecular dynamics by Iachello, Levine, Alhassid, Giirsey and their collaborators 
(see the book [5] for a survey of theory and applications) . Most of the applications of 
spectrum generating algebras concerned exactly solvable Hamiltonians whose spec- 
trum could be completely determined by algebraic means. An intermediate class of 
spectral problems are those for which only a finite part of the point spectrum can be 
calculated by algebraic methods, but possibly not the whole spectrum. The usual 
example of these type of problems is the sextic oscillator in which the potential is 
an even polynomial of degree six whose coefficients depend on a parameter n. For 
each positive integer value of n the Hamiltonian is shown to preserve an (n + 1)- 
dimensional subspace of Z^QR). It is clear from Sturm-Liouville theory that the 
Hamiltonian with the sextic potential has an infinite number of bound states, but 
only n + 1 of them belong to the so called algebraic sector. It was soon realized that 
the sextic oscillator was just the first example of a large class of systems having this 
property, and the works of Shifman, Turbiner, Ushveridze [6-8] and their collabo- 
rators initiated the study of the mathematical properties and physical applications 
of this new class of spectral problems, which they named quasi-exactly solvable. 
However, even in the simplest case of one spatial dimension, there is no way to 
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ascertain whether a given potential is quasi-exactly solvable or not, so two methods 
were proposed to construct large families of quasi-exactly solvable problems: the 
Lie-algebraic method [7] and the analytic method [9]. 

The idea underlying the Lie-algebraic method is to use results from the represen- 
tation theory of Lie algebras to ensure that a given Hamiltonian preserves a certain 
finite dimensional subspace of functions. In one dimension, the only algebra of first 
order differential operators with finite dimensional is sl(2), whose generators: 

Jn = z 2 D z — nz 7 

(1) J° n = 

Ja = D z 

preserve the (n + l)-dimensional linear space of polynomials in the z variable of 
degree less than or equal to n: 

(2) V n (z) = (l,Z,Z 2 ,...,Z n ). 

The most general second order differential operator that preserves V n (z) can be 
written as a quadratic combination of the generators JTJ of sl(2) 

(3) H = Y, C «* J n J l + Y, C " J n+ C *' 

where the indexes run over the values a,b — ±,0. The operator H is said to be 
Lie algebraic and sl(2) is often referred to as a hidden symmetry algebra. The 
operator H does not have in general the form of a Schrodinger operator, but it 
can be transformed into a Schrodinger operator by a change of variables and a 
conjugation by a non-vanishing function. Such a transformation always exists in 
one spatial dimension, but in more dimensions the equivalence problem remains 
open [10,11]. This has been a long standing obstacle to classify multidimensional 
quasi-exactly solvable Hamiltonians, where only a few families are known [12, 13]. 

One important fact lies at the core of the Lie-algebraic method: it is ensured 
by Burnside's classical theorem that every differential operator which leaves the 
space V n (z) invariant belongs to the enveloping algebra W(sl(2)), since V n is an 
irreducible module for the sl(2) action. This is probably the reason for the relative 
success of the Lie-algebraic constructions in the context of quasi-exact solvability, 
up to the point that the terms Lie-algebraic and quasi-exactly solvable are often 
used as synonims in the literature. 

However, there are many other finite dimensional polynomial spaces which are 
not irreducible modules for the sl(2) action, and in these cases there might be 
non-Lie algebraic differential operators which leave the space invariant. The first 
exploration of this type was done by Post and Turbiner [14], who classified all sec- 
ond order differential operators which preserve a polynomial space generated by 
monomials. In their work they did not use Lie algebras but other considerations 
based on the grading of the operators and basis elements. It was later realized that 
these differential operators can also be transformed into Schrodinger form thereby 
providing new examples of quasi-exactly solvable operators which are not Lie alge- 
braic. In this direct approach to quasi-exact solvability [15] more general polynomial 
spaces are considered and the set of differential operators that preserve them are 
investigated without any reference to Lie algebras. It was later shown that exactly 
solvable potentials exist which are not Lie algebraic, but can be obtained from a 
Lie algebraic potentials via a state-adding Darboux transformation [16]. Further 
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research along this line suggests to regard a Darboux or SUSY transformation not 
just as a transformation on the potential and eigenfunctions, but as an algebraic 
deformation of an infinite polynomial flag [17]. Non Lie-algebraic potentials appear 
also in the recent work of Gonzalez-Lopez and Tanaka in the context of supersym- 
metry [18]. 

In contrast with the Lie algebraic method, the analytic method has received less 
attention over the last decades. However recent results suggest that it is a more 
suitable approach to encompass this new type of quasi-exactly solvable systems 
which do not fit in the Lie algebraic scheme, [19]. 

In this paper we would like to illustrate the relevance of studying quasi-exact 
solvability beyond the Lie algebraic approach by providing two suitably chosen 
examples: in Section we show how more energy levels can be obtained from a 
Lie algebraic potential which cannot be obtained in the sl(2) approach. In Section 
we show an example of a potential which is quasi-exactly solvable but not Lie 
algebraic. 

2. Calculating more algebraic energy levels of a Lie algebraic 

potential 

In this Section we expose the first reason to study quasi-exactly solvable problems 
beyond the Lie-algebraic approach. We present a potential which is known to admit 
an sl(2) algebraization, and therefore allows for some of its energy levels to be 
calculated algebraically all belonging to the even sector. For this same potential, 
we show that a different algebraization exists in which the algebraic Hamiltonian 
cannot be expressed as an element of the enveloping algebra of sl(2). This new 
algebraization allows to calculate all the previous levels, plus some extra ones. 

Consider the following Schrodinger operator 

(4) H = —D xx + 2A 2 cosh(2x) + 4An cosh(x) - isech 2 (x/2), 

where A < is a real parameter and where n is a positive integer. This potential is 
known to be Lie-algebraic and it appears for instance in the classification performed 
by Gonzalez-Lopez, Kamran and Olver in [20]. The sl(2) algebraization of this 
potential is achieved by the transformation 

T = /i(x)^ 1 Hfi(x), 

with the following choice of gauge factor and change of coordinate 

(5) n(x) = e 2Acosh ^sech(x/2), 

(6) z = -sinh 2 (a;/2). 

Up to an additive constant, the transformed operator T becomes 

(7) T = z(l — z) D zz + (8Az(z -1) + -)D Z - 8Anz 

which is easily seen to preserve the space V n {z) defined in J2J). The explicit qua- 
dratic combination of T in terms of sl(2) generators is (again up to an additive 
constant) 

1 n 

(8) T = -( J°) 2 + - { J°, J-} + 8A J+ + (1-8A- n) J° n + - J- , 
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As a consequence of this sl(2) algebraization we can calculate dimV n — n + 1 
energy levels of the Hamiltonian, by diagonalizing the corresponding matrix of the 
restricted action of T to V n . The (n + 1) algebraic eigenfunctions have the form 

(9) ^2k(x) = fi(x) Pk (-smh 2 (x/2)), k = 0,...,n 

where pk is one of the (n + 1) polynomial eigenfunctions that T is ensured to have. 
All the eigenfunctions obtained via the Lie algebraic method correspond to the even 
sector. 

Up to this point, all these results are well known. However, the same Hamiltonian 
remarkably admits a different algebraization characterized by the following gauge 
factor and change of variables: 

(10) £(*) = e 2Acosh(:E) e-™ ;E sech(a;/2), 

(11) w = e x . 

The transformed operator T = [x{x)~ l H fi{x) becomes 

(12) f = -w 2 D ww +( -2Aw 2 + 2nw + 2(A - 1) + ——] D w +4Anw+-^—, 

\ 1 + w ) 1 + w 

where an additive constant has been dropped. It is straightforward to realise that 
due to the rational coefficients in the previous expression, the operator T docs 
not belong to the enveloping algebra of sl(2). Nevertheless it preserves a 2n- 
dimensional subspace generated by polynomials in the variable w that we shall 
denote as S^' ^ (w) and define as 

(13) St~ 1} = {?£ V 2n : p'(-l) + np(-l) = 0} 

This polynomial subspace is an example of an exceptional polynomial module, the 
name exceptional being due to the fact that the space of second order operators 
that leave it invariant have a rich structure [15,19]. An explicit basis of £^~ 1 \w) 
can be given in the following manner 

(14) 4T 1} M = (1 - n(w + 1), (w + l) 2 , (w + l) 3 , ...,(w + l) 2n ) 

This second algebraization allows us to compute 2n algebraic eigenfunctions of the 
Schrodinger operator which will be of the form 

(15) VfcO) = [i(x)p k (e x ), 

where pk(w) is one of the 2n polynomial eigenfunctions that T is ensured to have. 
We observe that these eigenfunctions can be both odd and even, as opposed to the 
sl(2) algebraic ones, which were only even. 

We thus obtain n+1 algebraic eigenfunctions via the sl(2) algebraization and 2n 
algebraic eigenfunctions via the exceptional module algebraization. The situation 
is depicted schematically in Figure ^ It is maybe worth to explore the relation 
between the two different algebraizations. It is known in the theory of Lie algebraic 
problems that under linear fractional transformations the Lie-algebraic character 
is preserved [20]. More specifically, a change of variable by a linear fractional 
transformation, 

z = a \ y 7 , aS-(3-f = l 7 
7z + d 
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sl(2) algebraization Exceptional 
algebraization 



Figure 1. The algebraic energy sector coming from both alge- 
braizations of potential with n = 4. 



together with an appropriate gauge transformation preserve the space V n since 
(16) p(z) G V n (z) \ * p(~z) = ( 7 5 + S) n p (^±4) G . 



jz + S 

Moreover, if T is Lie-algebraic and preserves V n (z), then 

(17) f = (-jz + ay n T(--fz + a) n 

is also Lie-algebraic and preserves V n {z). We see thus that projective transforma- 
tions do not change the sl(2) character, in fact they amount to a linear transfor- 
mation of the generators of sl(2) [20]. However, from expressions (|5|l.lfH|). (|10|) and 
(fllf) it follows that the relation between the two algebraic operators T and T in the 
previous example is given by 

(18) Z = 4^' 

(19) T = W - n fw n , 

which is not a projective transformation. 

We have performed an explicit computation to show a few of the algebraic eigen- 
functions of the Hamiltonian Q}. For illustrative purposes it suffices to fix the 
integer parameter n to some low value. For the values n = 3 and A = —0.4 the 
potential (@J has the shape of a double- well potential shown in Figure |21 

2.1. sl(2) algebraization. For n = 3 the operator T in Q) satisfies 

TV 3 (z)cV 3 (z) 

which allows to compute four algebraic eigenfunctions via the sl(2) algebraization. 
The action of T relative to the canonical basis of V 3 is given by the following 4x4 
matrix: 

-24A 
_| 1/2 -8A -16A 
3 -2 - 16A -8A 
15/2 -6-24A 



T 



V 3 
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For A = —0.4 we have calculated four algebraic eigenfunctions of (@J which have 
the form 



(20) 



fj,(x)p k (sinh -), 



fc = 0,1,2,3, 



where p k € Vziz) is one of the four eigen-polynomials of T and the gauge factor 
/i(x) is given by 0. These four eigenfunctions have been plotted in Figure along 
with their corresponding energies, which have been normalized so that the energy 
of the ground state is zero. 

2.2. Exceptional module algebraization. From the second algebraization we 
know that the Schrodinger operator H is conjugate via the change of variables 
and gauge factor i|10fl to the operator T which preserves an exceptional polynomial 
module. For n = 3 this space is spanned by 

(21) £i 3 '- 1} (w) = (-2 - 3w, (w + l) 2 , (w + l) 3 , (w + l) 4 , (w + l) 5 , (w + l) 6 ) , 
and the matrix of the action of T relative to this basis is 



12A 


-30A 














\ 


2 


10 - 4A 


8,4 














-6 


12 


6A 













-4 


-2 


12 + 4A 


4A 













-10 


6 


10 + 8A 


2A 













-18 


18 


6 + 12A 


J 
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V 

For the same value A = —0.4 we have calculated six algebraic eigenfunctions of the 
Hamiltonian |0J which have the form 

(22) *l> k (x)=fi(x)p k (e x ), k = 0,1,2,3,4,6. 

where p k £ f g ' is one of the six eigen-polynomials of T and the gauge factor jj,(x) 
is given by (|10() . These six algebraic eigenfunctions along with their energies are 
shown in Figure [3] The first thing to note is that we obtain all the even eigenfunc- 
tions from the sl(2) algebraization plus two extra odd eigenfunctions corresponding 
to the first and third excited states. However, the cigenfunction corresponding to 
the fifth excited state is not present in the algebraic sector. For arbitrary n it seems 
that there is always a gap in the algebraic sector just below the highest energy in 
the exceptional module algebraization, [15]. Although all evidences show that this 
is true in the general case, there is yet no proof of this result. 

3. A non-Lie algebraic quasi-exactly solvable potential 

In the previous section we have seen how the exceptional module algebraization 
can provide more energy levels from a Lie-algebraic potential. In this Section 
we will show that some Schrodinger operators only admit an exceptional module 
algebraization and not the traditional sl(2) one. We can therefore construct new 
quasi-exactly solvable potentials on the line, which are not in the classifications of 
QES potentials in [7,20], since those classifications deal only with potentials that 
admit an sl(2) algebraization. In this Section we show one first simple example 
of this phenomenon, postponing the full classification of these new potentials to a 
future publication, [19]. 

Consider the following Schrodinger operator: 



(23) 



H 



V(x) 
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with potential 

(24) V(x) = A 2 x 6 + 2ABx A + (B 2 + [{-If - 4n]A)x 2 + 4 



{x 2 + l) 2 ' 



If the last rational term were absent, this potential would be the well known quasi- 
exactly solvable sextic potential discussed for instance in [6] . The potential is always 
even, so its eigenfunctions will have well defined parity. However, as it happens with 
the sextic, only even or odd eigenfunctions but not both appear in the algebraic 
sector. This does not exclude in principle that other algebraizations exist in which 
both even and odd eigenfunctions are obtained. In fact, the choice of p = gives a 
potential with algebraic even eigenfunctions while p = 1 corresponds to potentials 
with odd algebraic eigenfunctions. For arbitrary values of A and B the Hamiltonian 
does not preserve any finite dimensional subspace, but for the following values: 



(25) A 



-iy 



n — a V 2 



(26) B = -^—(3a-2n+(-l) p (^-l) 

n — a \ 2a 

where a is an arbitrary real parameter, the Hamiltonian H admits an exceptional 
module algebraization. More specificallly, the transformation T = fi(x)~ 1 H fx(x) 
with 

(27) (j,(x) = -^_exp(-^x 4 -f:r 2 ), P ={0,1} 

x z + 1 4 2 

(28) w = x 2 

transforms H into the algebraic operator T given by 

T(w) = 4(J 4 + J 5 - A,h) 
up to an additive constant. The operators J4, J5 and Jq are given by 

(29) J 4 = (w-l)D ww + (a{w-l)-l)D w , 

(30) J 5 = D ww + 2(a —)d 



W — Is W — 1 

(31) J 6 - (w-l)(w-l--^D w -n(w-l), 

and every one of them leaves invariant the exceptional module S^'^iw) given by 

(32) 4 M) M = (1 - a(w - 1), (u, - l) 2 , (w - l) 3 , ...,(w- 1)") 

The sextic potential (|24|l has eigenfunctions in L2(M) provided that A < 0, or 
A = and B < 0. This last case corresponds in fact to an exactly solvable 
potential, which can be obtained by a SUSY transformation from the harmonic 
oscillator [16]. The algebraic operator preserves in this last case an infinite flag of 
exceptional polynomial subspaces. 

We show some explicit eigenfunctions of this potential corresponding to the even 
sector by setting n = 4 and p — 0. The matrix of the action of T with respect to 
the basis ^ of £ { 4 a,1) (w) is 
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T 





/ -15a 2 -a + 6 
a — 4 

-8 


6a 2 (l-2a) 
a-4 
-19a 2 +47a-10 




4a(2a— 1) 






2a(2a — 1) 


1) 




V 


a-4 

4(4a - 3) 
16 


a-4 
-23a 2 + 79a-18 




a-4 
8(3a - 4) 


a-4 
-27a 2 + llla- 
a-4 



For a = 0.8 the four polynomial cigenfunctions of T are approximately 



3.543u/ 



Po(w) = 

Pi(w) = 

P2(w) = 

p 4 (w) = 0.329 - 0.538u; + 0.114w 2 - 0.008u; 3 + 0.0001w 



0.431 

0.965- 

0.647- 



2.555w 
2.799w 
0.230w 



2.681u/ 
0.369w 2 



2.071w 3 
0.776w 3 
0.052w 3 



0.325w 4 , 
0.055w 4 , 
0.0017W 4 , 



where polynomial p^ has k zeros. The corresponding energies are: 

E a = 0, 
E 2 = 6.058, 
E 4 = 11.148, 
E s = 15.129. 

Therefore the four even eigenfunctions of l|24|) coming from the exceptional module 
algebraization are: 



(33) 



^2k{x) = H{x)p k {x 2 



k = 0,1,2,4. 



where the gauge factor ^i(x) is given by (|T7J) with p = 0. For this choice of param- 
eters a, p and n the potential, gauge factor and eigenfunctions have been plotted 
in Figure 0] We observe once again that we obtain eigenfunctions with zero, two, 
four and eight zeros, but not one with six zeros. For arbitrary n the eigenfunction 
with 2n — 2 zeros would be missing form the algebraic sector. The calculations 
for the odd case are similar: setting p = 1 in we have calculated n = 4 odd 
eigenfunctions which are shown in Figure [SJ 



4. Discussion 

By means of two examples we provide further arguments to motivate the study 
of quasi-exactly solvable problems beyond the Lie algebraic approach. Even in the 
simplest case of one dimensional quantum mechanical problems the classifications 
performed in [7, 20] do not cover all quasi-exactly solvable potentials. Since Lie 
algebraic potentials are only a subclass of the potentials with partial algebraization 
of their spectrum, it would be desirable to find a precise mathematical character- 
ization of the latter. This problem is related to the so called generalized Bochner 
problem [21], that of finding the most general differential operator that preserves a 
general finite dimensional space of polynomials. All the new quasi-exactly solvable 
potentials obtained recently are based on the exceptional modules, but other poly- 
nomial subspaces might exist which do not have a monomial basis and yet have a 
rich structure. In conclusion, the class of Hamiltonians with partial algebraization 
of their spectrum is larger than it was previously thought, and the new results call 
for new theoretical developments on this field. 
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Potential V(x) 



n{x) 




Gauge factor fi(x) 



E„=0 




1 2 3 X 



-3 -2 -1 

Ground state ipo(x) 



E,= 2.853 




Second excited state ip2(%) 



E 4 =7.183 




E„=1 3.740 




Fourth excited state 1/14(3;) 



Sixth excited state tpQ(x) 



Figure 2. The four algebraic eigenfunctions obtained through 
sl(2) algebraization of potential Q with n = 3 and A = —0.4 
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Figure 3. The six algebraic eigenfunctions obtained through the 
exceptional module algebraization of potential J3J with n = 3 and 
A = -0.4 
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Fourth excited state Eighth excited state tpe(x) 



Figure 4. Four even eigenfunctions of the modified sextic poten- 
tial H24f> with a — 0.8, n — 4 and p = 0, corresponding to the 
exceptional module algebraization. 
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First excited state ipi (x) Third excited state ips{x) 




Fifth excited state i/> 5 (x) Ninth cxcitcd state 

Figure 5. Four odd eigenfunctions of the modified sextic poten- 
tial H24fl with a — 0.8, n — 4 and p = 1, corresponding to the 
exceptional module algebraization. 



